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Motivated by a recent experimental report 1 claiming the likely observation of the Majorana mode 
in a semiconductor-superconductor hybrid structure^^, we study theoretically the dependence of 
the zero bias conductance peak associated with the zero-energy Majorana mode in the topological 
superconducting phase as a function of temperature, tunnel barrier potential, and a magnetic field 
tilted from the direction of the wire for realistic wires of finite lengths. We find that higher tem- 
peratures and tunnel barriers as well as a large magnetic field in the direction transverse to the 
wire length could very strongly suppress the zero-bias conductance peak as observed in Ref.[J|. We 
also show that a strong magnetic field along the wire could eventually lead to the splitting of the 
zero bias peak into a doublet with the doublet energy splitting oscillating as a function of increas- 
ing magnetic field. Our results based on the standard theory of topological superconductivity in a 
semiconductor hybrid structure in the presence of proximity-induced superconductivity, spin-orbit 
coupling, and Zeeman splitting show that the recently reported experimental data are generally con- 
sistent with the existing theory that led to the predictions for the existence of the Majorana modes 
in the semiconductor hybrid structures in spite of some apparent anomalies in the experimental 
observations at first sight. We also make several concrete new predictions for future observations 
regarding Majorana splitting in finite wires used in the experiments. 



I. BACKGROUND 



The search for solid state Majorana modes, which 
are localized quasiparticles with non-Abelian braiding 
statistics and a direct realization of the Majorana op- 
erator, has created a great deal of recent interest in 
the whole physics community^, partly because of the 
great mystique associated with Majorana and partly be- 
cause of the prospects for topological quantum compu- 
tation using solid state Majorana quasiparticles. The 
most experimentally promising among the many theo- 
retical proposals for solid state Majorana fermions is 
the semiconductor-superconductor hybrid structures^— 
where the ordinary s-wave superconducting proximity ef- 
fect induced in the semiconductor is modified by the in- 
trinsic spin-orbit coupling in the semiconductor and an 
externally applied Zeeman spin splitting. The solution of 
the resultant superconducting gap equations, often called 
the Bogulibov-De Gennes (BdG) equations, in the pres- 
ence of both spin-orbit coupling and spin splitting explic- 
itly shows a quantum phase transition in the nature of 
the superconducting phase as a function of the Zeeman 
splitting (which can be controlled by an external mag- 
netic field). The lower magnetic field part of the super- 
conducting phase is an ordinary non-topological s-wave 
superconductor with a suppressed gap due to the finite 
spin splitting which eventually gives way to a topologi- 
cal (effectively a p-wave) superconducting phase at large 
spin splitting when the original proximity-induced s-wave 
gap is completely suppressed. The Majorana mode ap- 
pears naturally as a stable zero-energy solution inside 
the gap of this topological superconducting phase. It has 



been shown that the existence of the Majorana quasi- 
particle would give rise to a conductance peak at zero 
energy inside the nominal p-wave superconducting gap 
of the topological phase. The experimental observation 
of this predicted zero energy conductance peak inside 
the superconducting gap at finite magnetic field in the 
semiconductor-superconductor hybrid structure is there- 
fore believed to be the necessary evidence in support of 
the solid state Majorana modei. 

It is useful to compare our semiconductor nanowire 
model for the Majorana modes with the original ID 
Majorana wire model introduced by Kitaev ("the Ki- 
taev model"). In the Kitaev model, the topological p x - 
wave superconductivity is assumed explicitly by making 
the system a spinless p-wave superconductor. In fact, 
Kitaev in his paper— goes on to discuss at quite some 
length the possibility of realizing such a topological ID 
spinless p-wave superconductivity in nature, concluding 
that the most suitable systems for realizing such spin- 
less ID p-wave superconductivity are quasi-lD organic 
superconductors involving charge and spin density wave 
ground states, which were already proposed as possible 
Majorana-carrying systems by Sengupta et ali^ who in 
fact also pointed out that the zero bias conductance peak 
associated with such Majorana modes would be quan- 
tized. It is a curious historical coincidence that Kitaev 
in his original paper explicitly rejected both spin-orbit 
coupling and Zeeman splitting as possible physical mech- 
anisms for producing ID spinless p-wave superconduc- 
tivity! The current activity in ID semiconductor Ma- 
jorana wires started from the work of Sau et ali^ who 
showed that in 2D semiconductors with proximity s-wave 
superconductivity induced by a nearby ordinary s-wave 
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superconducting metal, the topological superconductiv- 
ity (" 2D chiral p-wave superconductor" ) would naturally 
arise in the presence of spin splitting and spin-orbit cou- 
pling provided that the spin-splitting is large enough to 
overcome the trivial s-wave superconductivity. It was 
soon realized^— that this 2D chiral topological super- 
conductivity can easily be modified to ID helical topo- 
logical superconductivity by considering ID semiconduc- 
tor nanowires rather than 2D semiconductor heterostruc- 
tures with the spin splitting being introduced by an ex- 
ternal magnetic field. This helical ID p-wave supercon- 
ductor is an effective realization of the ID Kitaev model, 
but it is not identical to the Kitaev model. For example, 
the topological superconductivity in the ID nanowires 
exists only above a finite value of the spin splitting with 
only trivial non-topological superconductivity existing at 
weaker magnetic field values. Also, the superconducting 
gap depends explicitly on the spin splitting and the spin- 
orbit coupling in the semiconductor nanowires instead of 
being just effective theoretical parameters. 



II. INTRODUCTION 



In a recent presentation^, the likely experimental ob- 
servation of the theoretically-predicted^— zero-energy 
Majorana modes in semiconductor nanowires, in close 
proximity to an ordinary (i.e. s-wave) superconductor 
and in the presence of an external magnetic field ap- 
plied along the wire, has been reported. This experi- 
mental observation of the predicted zero bias peak and 
the associated implication that this may finally be the 
real-life evidence for the existence of the elusive (and 
so-far purely theoretical) Majorana mode have created 
tremendous excitement in the general scientific commu- 
nity. Given the great significance of the possible exper- 
imental discovery of the Majorana mode, it is therefore 
of utmost importance to ensure that all aspects of the 
experimental discovery in Ref.[l] are indeed consistent 
with the theoretical expectations and there are no loose 
ends. This is particularly true in view of the facts that 
the experimental observation precisely follows theoretical 
predictions for the existence of the Majorana zero energy 
mode in semiconductor-superconductor hybrid structures 
and that zero bias peaks could arise in superconductors 
and semiconductors from a variety of physical effects such 
as Andreev or Shiba bound states, Kondo resonances, 
etc. Our work presented in this paper aims at a qualita- 
tive understanding of several interesting features of the 
experimental observation to ensure that the observation 
is consistent with the expectations of the BdG theory 
for the Majorana mode in semiconductor-superconductor 
structures. 

In the current work, we theoretically investigate some 
of the peculiar aspects of the experimental observa- 
tions in Ref. [l[ which were not directly or explicitly pre- 
dicted earlier— _ —i2r— in the extensive theoretical work 



leading to the experimental observation of the Majo- 
rana mode although some of the aspects we discuss in 
our work were implicit in the theory. The experimen- 
tal observation specifically concentrates on the study of 
a zero-bias-conductance peak (ZBCP) in the current (I)- 
voltage(V) differential conductance (^) of the tunnel- 
ing spectroscopy of an InSb nanowire on superconducting 
NbTiN, which manifests itself only in the presence of an 
external magnetic field B x (> 0.1 T) oriented along the 
wire (taken to be the x-axis in this paper) . The existence 
of this ZBCP in Ref.[l[ for B x ^ has been claimed to be 
the verification of the theoretical prediction for the exis- 
tence of the zero-energy Majorana mode^— in the wire. 
The zero-energy Majorana mode exists as localized quasi- 
particles at the ends of the superconducting wire and is a 
manifestation of the system being in a chiral p-wave topo- 
logical superconducting (TS) phase as envisioned more 
than a decade ago^^. The theory predicts^— the 
presence of the TS phase for V x > V c = -\/A 2 + /i 2 , 
i.e. B x > B c = V c /giiB with V x = g/iBB x being the 
Zeeman field in the wire associated with B x , and A,// 
are the superconducting gap and the chemical potential 
respectively in the wire. For B x < B c (or V x < V c ) 
the system is in an ordinary non-topological (i.e. s- 
wave) superconducting phase (NTS) which in the pres- 
ence of the finite Zeeman splitting V x makes a topologi- 
cal quantum phase transition^— Zr— to the TS phase for 
V x > V c (i.e. B x > B c ). The TS phase has the Ma- 
jorana modes localized at the ends of the wire and the 
associated ZBCP at zero energy in the middle of the su- 
perconducting gap. The NTS phase on the other hand 
has no structure, except perhaps some Andreev bound 
states (ABS) at generic non-zero energies, within the su- 
perconducting gap. The existence of a robust ZBCP in 
the differential tunneling conductance has therefore been 
predicted^— &2 r 11 1 13 to be the necessary condition for the 
observation of the Majorana mode, and its observation in 
Ref.[l[ is an important experimental milestone providing 
perhaps the first definitive signature for the Majorana 
fermion in a solid state system. 

Given the key importance of the subject matter, 
namely, the possible experimental discovery of the emer- 
gent Majorana mode in the topological superconductor 
system, it is somewhat disconcerting that some of the ob- 
served experimental features are unexpected and some- 
what inconsistent with the existing theoretical predic- 
tions in the literature although most of the findings in 
Ref.pJ are, in fact, completely consistent with the the- 
oretical predictions (e.g. the existence of ZBCP only 
above a critical value of B x ). We concentrate in the cur- 
rent work on three features of the experiment which, in 
our opinion, require special attention: (1) the observed 
ZBCP is much (by more than an order of magnitude) 
weaker than the predicted^ ' 1 1 1 13 ~ — canonical quantized 
value of 2e 2 /h expected for the Majorana zero energy 
mode; (2) a peculiar and unexpected splitting of the 
ZBCP at high values of V x (for B x > 0.5 T) observed 
in Ref.[l]; (3) the predicted behavior of the ZBCP in 
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the presence of an additional transverse component V y 
of the Zeeman field associated with an applied magnetic 
field component B y (= Vy/g/Js) along the direction of the 
spin-orbit coupling field (y-axis) which is known 16 to be 
transverse to the length to the wire. Of the three issues 
theoretically considered in this work, the first two are di- 
rectly motivated by the experimental data presented in 
Ref.[l| where a strongly suppressed ZBCP (with a differ- 
ential conductance value substantially below 2e 2 / h ) and 
a splitting of the ZBCP into a doublet at high values of 
B x are both observed. Item (3) in our work is alluded to 
in Ref.JIj], and our work here provides the numerical re- 
sults for the expected experimental observation when the 
applied in-plane B field is tilted at an angle 9 to the wire 
length direction, i.e. (B x ,B y ) = (B cos 9, B sin 9) , which 
gives (V x ,V y ) = (g/iBB cos 9, g[iBB sin 9) where g, [Ib are 
the Lande g-factor and the Bohr magneton respectively. 

The most important new qualitative feature of our the- 
oretical work presented in this work is considering realis- 
tic experimental systems, in particular, finite wire lengths 
and finite temperatures as well as finite tunneling barrier 
heights. We find that the finiteness of the nanowires is 
a fundamental constraint in the ideal realization of the 
Majorana zero energy mode, and most of our interesting 
and important results arise directly from our keeping wire 
lengths finite as in the experimental systems. The reason 
for the qualitative importance of the finite wire length is 
rather obvious. The external magnetic field suppresses 
the superconducting gap, both in the NTS and the TS 
phase, thus enhancing the coherence length which varies 
inversely as the gap energy. When the enhanced coher- 
ence length becomes comparable to the wire length, the 
two end Majorana modes start 'seeing' each other, lead- 
ing to an energy splitting. This effect turns out to be of 
qualitative importance, affecting for example the conduc- 
tance quantization of the Majorana ZBCP even at T = 
in contrast to the ideal infinite wire case where the ZBCP 
is always quantized at 2e 2 /h in the T = limit. 

The quantization of the ZBCP predicted for Majorana 
fermions^ i 1:L i 13 ' 15 is a result that is valid only in the zero- 
temperature limit. In contrast, the high-temperature 
limit has a more conventional resonant scattering Fano- 
form — i 15 ' 17 with a height proportional to Y /ksT and 
a width proportional to the thermal energy ksT, for 
ksT 3> r with T/h being the tunneling rate between 
the Majorana bound state and the lead. The tunneling 
rate T/h depends on the transparency of the tunneling 
barrier, and is not easily experimentally controllable- in 
fact, the barrier transparency and hence the tunneling 
rate is simply unknown in the experimental situation. 
In the part of this work which focuses on Item (I), we 
show how the tunneling conductance crosses over from 
the low-temperature to the high-temperature limit and 
establish that for reasonable parameters, it is indeed pos- 
sible to have a dramatic suppression of the ZBCP as seen 
in Ref. 1]. In addition, from our numerical calculations 
for the realistic parameters, we have found that even at 
T = 0, the ZBCP can be suppressed below its quan- 



tized value for sufficiently small tunneling rate T because 
of finite size effects. In particular, such a T — sup- 
pression of ZBCP happens when T becomes compara- 
ble to the splitting between the end Majorana fermions, 
which may be the case for the few micron long wires in 
Ref.[l| together with the tunneling rates T inferred from 
the measured conductance. The splitting between end 
Majoranas is invariably present in real wires where the 
two end Majoranas have some finite overlap, leading to 
a lifting of their precise zero energy status. The splitting 
of the ZBCP as a function of V x , which is discussed as 
part of Item (2) in this paper is a finite size effect, which 
is likely to be relevant for the experiments in Ref. [ij , but 
as far as we are aware has not been discussed in the 
literature except in the idealized situation^. The pre- 
dictions in the literature, which are restricted to infinite 
wires, show that the Majorana fermion must be robust 
for large Zeeman fields in the case of narrow wires where 
the inter-sub band spacing is much larger than the Zee- 
man splitting. The numerical results presented in this 
paper show that for finite wires, even in the narrow wire 
limit, the ZBCP is split for large V x . The splitting of 
the ZBCP arises from overlap of the Majorana fermion 
wave-functions as has been previously discussed in the 
context of p-wave superconductors^. To the best of our 
knowledge, our work is the only existing result in the lit- 
erature for the Majorana splitting in finite wires in the 
presence of spin-orbit coupling and Zeeman splitting. Fi- 
nally in the context of Item (3) , we discuss the effect of 
the angle of the Zeeman potential on the ZBCP. Consis- 
tent with previous theoretical work 3 , which shows that 
the proximity-induced quasiparticle gap vanishes in the 
wire for V y > A, we find that the ZBCP vanishes above 
a threshold value of V y . 

We emphasize that our goal in this work is not to de- 
velop a new theory for the Majorana modes in the semi- 
conductor nanowires, but to extend and generalize the 
existing theory^— leading to the prediction of the Majo- 
rana mode in semiconductor hybrid structures, to finite 
wire lengths, temperatures, and barrier heights as well as 
to large Zeeman splitting (along arbitrary directions also) 
to see whether the features arising out of the standard 
theory are consistent with the observations of Ref.[l[. 
Of course, theory by itself cannot answer the question 
whether the observations in Ref. [I[ are really the isolated 
Majorana modes predicted by theory^— or are some dis- 
order or multi-mode wire effect. Instead our goal is to 
show that some of the qualitative features of the experi- 
ment, which are not immediately consistent with the re- 
sults explicitly predicted in previous works^— , are indeed 
consistent with the experiment, therefore, it is impor- 
tant to discern which features of the experimental data 
are in direct agreement with the theory extended to in- 
clude finite wire length, finite temperature, finite barrier 
height, and large Zeeman field. Such a comparison, as 
carried out in this work, between the extended standard 
theory and the data of Ref. [I[ will help indicate whether 
the experiment really is consistent with a Majorana in- 
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terpretation and also whether new ideas are essential to 
understand some features of the data. As discussed in 
more detail in Sec. VI, we intentionally avoid discussing 
issues such as disorder and the multi-band effects, since 
these effects would complicate the conclusions by intro- 
ducing more unknowns into the theory. Moreover, as 
elaborated further in Sec. VI, the experimental data in 
Ref.[l[ appears to use wires with relatively long mean- 
free paths (300 nm) and have large sub-band spacings so 
that such effects are unlikely to change our conclusions in 
a qualitative way. In this context, it is also important to 
emphasize that the finite wire length automatically im- 
plies a crossover in the behavior of the Majorana mode 
when the magnetic field induced suppression of the su- 
perconducting induced gap leads to the coherence length 
becoming comparable to or larger than the wire length. 
This physics is qualitatively new in our work (although 
it is implicit in the earlier works) because a finite wire 
length allows the two Majorana modes to overlap with 
each other leading to a splitting which is by definition im- 
possible in an infinite wire at any magnetic field. Since 
the experiments are always done in finite wires, our work 
provides a crucial extension of the standard theory in or- 
der to understand or interpret the experimental results 
even at a qualitative level. 



III. THEORY 



The physical system^ for studying Majorana 
fermions includes a strongly spin-orbit coupled semicon- 
ductor (SM), proximity-coupled to an s-wave supercon- 
ductor (SC) and imposed to a Zeeman field. Without 
loss of generality, we consider a finite ID SM nanowire 
along the x direction, the spin-orbit interaction an, be- 
ing along the y axis, and a Zeeman field V = {V x ,V y ). 
Also the wire is in contact with a superconductor, with 
proximity induced pairing strength A. The continuous 
BdG Hamiltonian for the system is 
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t z + V X (T X + V y a y +ia R d x a y t z + At x 



(1) 

\x is the chemical potential. The Pauli matrices er, r op- 
erate in spin and particle-hole space, respectively. For 
numerical convenience, it is standard to study the BdG 
equation in a discrete lattice tight-binding approxima- 
tion with no loss of generality and with fewer unknown 
parameters. Under the lattice approximation, we can 



map ([J) to a tight-binding model 

H = Y1 *y' c la c j<r - VicLct* (2) 

i,j,cr i,cr 
i,crcr' 

+ 4* ( V =° (i) Ox + V y (i) Oy) aa , c l(jl 

i 

The first contribution describes hopping or the kinetic 
energy along the length (i.e. the x-axis) of the wire, 
the second term represents the Rashba spin-orbit interac- 
tion, the third term is the Zeeman field and the last term 



shows the proximity-induced pairing term. cJ 



de- 



note electron creation and annihilation operators, respec- 
tively. We include only nearest-neighbor hopping with 
tij = —to and also include an on-site contribute ta = 2to 
that shifts the bottom of the energy spectrum to zero 
energy. The chemical potential /i is calculated from the 
bottom of the band. In the long wavelength limit, the 
tight binding model reduces to the continuum Hamilto- 
nian (TT]) with to — ti 2 /2m*a 2 and Rashba spin-orbit cou- 
pling an = aa with lattice constant a. In the numerical 
calculations we use a set of parameters consistent with 
the properties of InSb, as in Ref. [1], and choose the ef- 
fective mass to* = 0.015m e , spin-orbit coupling cur = 0.2 
meVA, and g^BjnSb = 1-5 meV/T. We also use A = 0.5 
meV and the length of the wire to be 4.5 fj,m. These pa- 
rameters are roughly consistent with the experimental 
conditions of Ref.[l[ although we are not interested in 
fine-tuning parameters for quantitative agreement with 
the data because there are far too many unknown pa- 
rameters in the experiment. We choose the tight binding 
numerical lattice parameter a = 15 nm, which is chosen 
so that the band- width 2to ^> V, A, /j, and thus the tight 
binding approximation itself would not introduce any ar- 
tifacts into our results. The length of the wire L — 4.5 /ito 
then corresponds to TV = 300 sites. We mention here that 
the SC proximity effect has now been observed by sev- 
eral groups in the SC/SM hybrid systems including both 
InSb nanowire o 1 ' 19 ' 20 and InAs nanowire a 21 ' 22 . Thus, 
the effective model given by EqfSl which is the starting 
point for our theory, is an appropriate model for the SM 
nanowire in the presence of spin-orbit coupling, Zeeman 
splitting, and proximity-induced superconductivity. 

To calculate the differential conductance measured 
from tunneling into the end of the superconducting 
nanowire, we have to study the current flowing into the 
wire contacted through a barrier region at one end with 
a lead by using the Blonders-Tinkham-Klawijk (BTK) 
formalism^ The main idea is to get the reflection and 
transmission coefficients by solving a BdG eigen-equation 
with some initial conditions. The current and conduc- 
tance can be expressed in terms of these coefficients. 
Given that the experiment involves many unknown pa- 
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rameters controlling many interfaces (e.g. SC/SM, tun- 
neling barrier/SM, various gates to control the barriers 
and chemical potential), our goal is to utilize the sim- 
plest possible model with the fewest number of parame- 
ters which would be capable of capturing the underlying 
qualitative physics of the experiment in finite wires. The 
BdG-BTK formalism provides the simplest qualitative 
basis for the theoretical modelling of the experimental 
system in Ref. [1]. 

More precisely, we start with the tight-binding Hamil- 
tonian @ with open boundary conditions so that the 
part i = 0, 1, ...mi are in the lead with no superconduc- 
tivity and the sites i = ttll + 1,...,N are in the super- 
conducting nanowire. The barrier region is modeled as 
a variation in the local chemical potential (M — > fM — U 
for the sites i — mi, — 2, ... , mi + 2, where U is the tun- 
nel barrier height. To calculate the reflection and trans- 
mission coefficients at energy E, we note that since the 
lead is normal (without spin-orbit coupling), the incom- 
ing mode can be taken to be purely electron-like. The 
reflected amplitudes can be normal r^.a^a' or anomalous 
fA.a,a'- Here a is the spin of incoming electron and a 1 
is the spin of the reflected electron and hole. The BdG 
equations that determine the reflected coefficients are 

{H m ,n - ES m . n ) = for m = 1, ..N 

n=0,N 



*n=0 = 




/ 





J,k h a 



The first equation is the BdG equation in the supercon- 
ducting wire, while the latter two equations express the 
wave-functions in the lead in terms of the reflection co- 
efficients. These equations need to be solved for both 
a =ti 4- • I n the above k e is the wave-vector in the lead 
at energy E so that 

k 2 
2m 

where fiiead is the chemical potential of the lead, while 
kh satisfies 

K h _ „ — _P 

2m 

The voltage-bias V of the lead determines the occupancy 
of the incident electrons. Electrons are incident on the 
superconductor from energy E = —fJ.i ea d to E = V. The 



states with normal reflection do not contribute to the 
current. Thus, the current will be 

7 = E r dEj2\rA,«,<r'(E)\ 2 

which implies a conductance (in unit of 2e 2 /h) 



dl_ 

dV 



= 5>,W (E = V)\ 2 . 



Technically for negative bias voltages V < 0, one needs 
to consider holes below the fermi-energy incident from 
the right with energy E — — V > and then Andreev re- 
flected at the interface and becoming electrons. However, 
such processes are related by particle-hole symmetry to 
Andreev reflection of negative energy electrons. There- 
fore, one does not need to calculate the Andreev reflection 
process separately since it is automatically incorporated 
in our BdG-BTK formalism. 

To generalize to finite temperatures all one needs to do 
is to broaden the conductance with the derivative of the 
Fermi-function in the usual manner 



G(V,T)= / deG (e) 



1 



4Tcosli ((!/-£) /2T) 



with Go (e) the zero temperature conductance at energy 
e. 

We present our numerical results for the calculated 
tunneling conductance as a function of bias voltage and 
other relevant experimental tuning parameters in Figs.l- 
5 to be compared with the experimental data of Ref.[l]. 
We first mention the fact that the (Majorana) properties 
of the real experimental system depend on (at least) ten 
independent parameters (many of which are unknown). 
These parameters include temperature, tunneling barrier 
(which in all likelihood depends on several unknown pa- 
rameters determined by the details of the interface and 
the control gates), Zeeman fields (V x and V y ), spin-orbit 
coupling, chemical potential, the induced superconduct- 
ing gap (which in turn depends on several parameters 
such as the semiconductor-superconductor hopping am- 
plitude, disorder, and the parent gap in NbTiN), the pa- 
rameters defining the ID confinement in the wire (which 
requires at least four independent parameters for con- 
finement along y and z directions), wire length (L) along 
the wire, and disorder (which by itself would necessi- 
tate several independent parameters for its description 
since in principle there could be long-ranged and short- 
ranged random impurities in the wire as well as interfa- 
cial roughness at the semiconductor-superconductor in- 
terface). No meaningful theory, beyond mere data fitting, 
can of course attempt to include all these parameters. 
In the current work we are interested in the fundamen- 
tal question of whether a minimal theoretical model can 
capture the basic qualitative findings of Ref.jlJ, and as 
such we ignore all the complications, concentrating on 
the single subband ID model in the absence of disorder 
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FIG. 1: (Two Left Panels) The differential conductance for a 
fixed Zeeman potential V = (1,0) meV and different tunnel- 
ing barrier U — 38 and 42 meV. The green dotted, red dashed 
and orange lines denote different temperature 0, 60, 120 mK, 
respectively. (Two Right Panels) Zero bias conductance peak 
as a function of temperature for Zeeman potential V = (1,0) 
meV and different tunneling barrier U = 38 and 42 meV, 
respectively. The peak decreases monotonously with tem- 
perature. (L = 4.5/im, spin-orbit coupling a = 0.2 meVA, 
induced pairing potential A = 0.5 meV) 



within a tight-binding BdG-BTK formalism. Our work 
should be construed as a zeroth order effective model of 
the experiment which should be the starting point for 
future quantitative and realistic theoretical studies when 
the details of the experimental parameters become avail- 
able. 



IV. RESULTS 



In FigQ] we show our calculated differential conduc- 
tance dl/dV as a function of the tunneling bias voltage 
V (which should not be confused with the Zeeman fields 
V x ,V y ) for two different tunnel barriers and three different 
temperatures for V x = 1 meV and V y = 0. (This choice 
of V = (V x , V y ) guarantees that the system is in the TS 
phase satisfying V x > \/A 2 + /i 2 for our choice of system 
parameters corresponding to Ref. [1].) In the third panel 
of FigfT] we depict 4^ of the ZBCP itself as a function 
of T for two values of U. These results manifestly estab- 
lish that the canonical quantized value of 2e 2 /h is clearly 
an unphysical theoretical limit achievable only as T — >• 
(and for low values of U). For reasonable values of U 
and T, our calculated value of ZBCP in FigJT] could eas- 
ily be one to two orders of magnitude smaller than 2e 2 /h, 
thus providing a satisfactory probable explanation for the 
weak strength of the ZBCP observed in Ref. [l| . We note 
that in our model, low values of U correspond to the more 
transparent barrier with higher conductance values- our 
numerical results clearly indicate that the ideal Majorana 



conductance of 2e 2 /h is unphysical even in the extreme 
low-temperature limit since the tunnel barrier U is never 
likely to be very small (or equivalently the transparency 
to the Majorana is not likely to be large). Our work leads 
immediately to the somewhat disappointing conclusion 
that the full Majorana spectral strength may never be 
achievable in the standard experimental set ups for mea- 
suring dl/dV since it would be very hard to get to the 
ideal limit of both T, U being very small (compared with 
the induced topological gap). Without our numerical re- 
sults, this would not be apparent. Since the conductance 
calculated for the larger barrier height in Fig.[JJis found 
to be significantly below the quantized conductance even 
at T = 0, it is clear that the suppression of conduc- 
tance in Fig. [JJis not purely the previously studied finite 
temperature effec t 11 ' 17 . Of course, as mentioned in the 
Introduction, this finding of ours, namely that the Majo- 
rana tunneling differential conductance is not quantized 
even at T = for finite U is valid only for finite wires 
since it is well-established that in the thermodynamic 
limit at T = the Majorana conductance is necessarily 
quantized to a value of 2e 2 /h. The key here is that in 
finite wires the two end Majoranas always overlap, and 
one must have the tunneling barrier low enough for this 
overlap to be unimportant in determining the conduc- 
tance quantization- thus for very long wires the conduc- 
tance quantization will be valid up to very large values 
of U, but for short wires, the value of U must be quite 
low. Our finding of the strong suppression of the ZBCP 
in short wires even in the T = limit due to the tunnel 
barrier effect (coupled with the Majorana overlap) is a 
new and somewhat unexpected result, which is consis- 
tent with the experimental finding of a saturation of the 
ZBCP at low temperatures^. While splitting of the MFs 
resulting from the finite length of the wire is expected 
to suppress the ZBCP because splitting of the MFs re- 
moves the zero-energy state responsible for the ZBCP, 
one would expect such splitting to result in a split ZBCP 
rather than one that is simply reduced in magnitude. A 
split ZBCP should be discernible in experiments at low 
temperatures and in this case the splitting of the MFs 
would be easy to determine from the experimental re- 
sults. Similarly, a reduction of the ZBCP from finite tem- 
perature would lead to a temperature dependenct ZBCP, 
whose width is comparable to fcgT. Our results plotted 
in Fig. [TJ show that for an intermediate regime of the 
barrier height, U, the finite overlap of the MFs can sup- 
press the ZBCP without leading to an observable split- 
ting of the ZBCP. Extremely large barrier heights U lead 
to peak width, which is dominated by temperature, and a 
temperature-dependent ZBCP height. Very low temper- 
atures T and high barrier heights could in principle lead 
to a regime where a split conductance is observed. On 
the other hand, low temperatures and a barrier with a 
transparency smaller than the finite size splitting would 
result in a ZBCP, which is at the same time indepen- 
dent of temperature and also have a ZBCP height that is 
significantly smaller than 2e 2 /h. Our results show that 
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FIG. 2: The differential conductance for a fixed V y — meV 
and different 14 = 0,0.5,1,2,3,4 meV. The red solid and 
blue dashed lines denote temperatures T — and 60 mK, re- 
spectively. At V x = 1 meV, the system is in the topological 
phase with quantized ZBCP. For length L = 4.5/im, the Zee- 
man splitting 14 > 1 meV reduces the gap, and leads to a 
stronger overlap between the MFs that results in a splitting 
of the ZBCP peaks. Additionally the finite temperature will 
decrease ZBCP. (Parameters used in the plot: U — 38 meV, 
L = 4.5/im, spin-orbit coupling a = 0.2 meVA, induced pair- 
ing potential A = 0.5 meV) 



finite size effects can suppress a ZBCP significantly below 
2e 2 //i, without the peak being temperature-independent 
or split. 

In Figl^l we show our calculated magnetic field or Zee- 
man splitting V x dependence (V y — 0) of dl/dV for two 
fixed temperatures keeping all other parameters fixed. 
The interesting result here in partial agreement with 
Ref.[l| is the splitting of the ZBCP for large (V x > 2 
meV) values of the Zeeman splitting. This ZBCP split- 
ting arises from the wire length (L — 4.5 /xm) being finite 
(both in our simulations and in Ref.[l|), which leads to 
the possibility of the two Majorana modes localized at 
the two ends of the wire to hybridiz o 18 i 24 causing the 
splitting of the ZBCP. (The double-peak structure of the 
ZBCP is apparent in Ref.[l] for the applied magnetic 
field above ~ 0.48 T where two peaks around zero bias 
can be seen in the data.) The Majorana hybridization 
effect (and consequently the splitting of the ZBCP) is 
exponentially suppressed for smaller values of V x (> V c ) 
still within the TS phase since the superconducting gap 
is large. With increasing V x , the gap is eventually sup- 
pressed as V~ x , which increases the coherence length, 
leading to an effective overlap between the Majorana 
modes localized at the two ends of the wire. There is 



nothing mysterious about the splitting of the ZBCP at 
high values of V x where the induced gap is small in the 
TS phase; this is expected- the important point is our 
finding that this happens in the same range of V x values 
in our theoretical modeling as it does in the experiment. 
Interestingly for the parameters of the problem in Fig 121 
we find the Zeeman induced splitting of the ZBCP to 
be only weakly dependent on Zeeman field. In partic- 
ular, consistent with the data in Ref.[l] our calculated 
ZBCP splitting (~ 0.02 meV) in Fig. 2 is much smaller 
than the applied Zeeman field (V x ^2 — 4 meV) caus- 
ing this splitting- this implies that the ZBCP splitting 
is not a trivial spin splitting either in our theory or in 
the experiment. As already mentioned above, the ZBCP 
splitting in the theory has its origin in the splitting of 
the Majorana zero energy mode due to the finite overlap 
of the two end Majorana localized wavefunctions over- 
lapping due to the finite length and the high field situ- 
ation, as predicted originally in Ref. [18(. We emphasize 
that the ZBCP splitting depends on increasing V x in a 
nontrivial manner and is not simply a Zeeman splitting 
going as linear in B, which serves to distinguish the high- 
field Majorana splitting from any run-of-the-mill Zeeman 
splitting arising in the experimental situation. 

The results presented in Fig. 2 clearly indicate that 
the ZBCP splitting (i.e. the Majorana splitting in finite 
wires) as a function of the Zeeman splitting is a non- 
monotonic function of the Zeeman field Vx, as was first 
predicted theoretically in Refs.[il|and [24j in a slightly 
different context. In fact, the splitting oscillates with the 
magnetic field V x since it is an oscillatory function of the 
superconducting coherence lengt h 18 ' . This feature was 
also observed in the numerical work of Ref. [l4[ ■ We dis- 
cuss this Majorana splitting in more details in discussing 
Fig. 6 below, but the oscillatory behavior is already ap- 
parent in the results for the various V x values shown in 
Fig. 2 for smaller values of V x . 

The observed high-field (> 0.48 T) splitting of the 
ZBCP in Ref.pJ could thus probably arise from a finite 
wire length effect in the high field regime. Of course a 
quantitative comparison with experiment of the exact na- 
ture of the splitting of the ZBCP would require a system- 
atic determination of the disorder and field configuration 
of the experiment. Also, the experimental data in Ref. [1] 
unfortunately end around 0.5 T value of the magnetic 
field, and it is therefore not possible right now to quan- 
titatively compare the field dependence of our numerical 
finding of the ZBCP splitting arising from the Majorana 
overlap with the experimental observation. Hopefully, 
our theory will motivate further high- field studies so that 
the dependence of the ZBCP splitting on V x can be com- 
pared between experiment and theory. 

It is know n 18 ' 24 that the Majorana splitting due to 
the overlap of two Majorana modes decays exponen- 
tially with their separation with the characteristic decay 
length being the superconducting coherence length. The 
main physics of the enhanced splitting with increasing V x 
found in Fig. 2 is thus the physics of the enhancement of 
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FIG. 3: The differential conductance for a fixed V x = 1 meV 
and different V y = 0.1,0.4,0.6 and 0.8 meV. The red solid 
and blue dashed denote temperatures T — and 60 mK, 
respectively. At V x = 1 meV, the system is in the topological 
phase with quantized ZBCP. V y will reduce the quasi-particle 
gap and hence suppress the ZBCP. Finite temperature also 
reduces the ZBCP. (U = 38 meV, L = 4.5/im, spin-orbit 
coupling a — 0.2 meVA, induced pairing potential A = 0.5 
meV) 



the coherence length due to the gap suppression by the 
external magnetic field. Such a splitting with increasing 
V x is thus expected at high field values, but the impor- 
tant thing to note is that the splitting is consistent with 
the experimentally observed ZBCP splitting in the data 
of Ref. [1], lending credence to the claim that Ref.[l] is 
really exploring Majorana physics. Note that in an infi- 
nite system, such a splitting would never arise and the 
ZBCP will never split no matter how large V x is. It is 
also known that in addition to the decay of the MF split- 
ting in finite systems, characterized by the separation 
between the MFs normalized to the superconducting co- 
herence length, the splitting also shows a characteristic 
oscillatory behavior as a function of the coherence length 
itself. In Fig. 2, V x is kept fixed in each panel, so these 
oscillations arc not apparent, but we discuss these MF 
splitting oscillations as a function of V x (through the en- 
hancement of the coherence length by V x ) in the Sec. V 
of the paper (see, e.g., Fig. 6(d) which explicitly shows 
the oscillation as a function of V x ) . 

Having provided reasonably realistic probable expla- 
nations for the two observed puzzling features of Ref. [l| , 
namely, the suppressed values of ZBCP (FigO} and the 
splitting of ZBCP (FigJ5J), we now consider the effect of a 
Zeeman field V y in the spin-orbit coupling direction (i.e. 
transverse to the wire length). If V y 3> V x , we expect 
the ZBCP to disappear even if the system remains in 
the TS phase in accordance with the invariant Pfaffian 
calculation 25 . This is because a large V y > A is known 
to suppress the quasiparticle gap rendering any end state 
completely delocalized across the wire4 In FigsJ3] and 0] 
we present our predicted results for the tunnel conduc- 
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FIG. 4: (Left Panel) ZBCP as a function of V y for a fixed 
V x = 1 meV at zero temperature denoted by the blue solid 
line. ZBCP shows a plateau for small extent of V y and is 
then suppressed by larger V y . The red dashed line denotes 
the quasiparticle energy gap. The peaks appear after gap 
closure. (Right Panel) Blue solid and red dashed lines denote 
different temperature T = 60 and 120 mK, respectively. At 
finite temperature the ZBCP is suppressed. [L = 4.5 iim, 
spin-orbit coupling a — 0.2 meVA, induced pairing potential 
A = 0.5 meV) 




FIG. 5: ZBCP as a function of angle 9 = tan' 1 (V y /V x ) for 
a fixed |V| = 1 meV at T = 0,60,120 mK denoted by the 
red dashed, blue dotted and orange solid lines, respectively. 
(L — 4.5 /jm, spin-orbit coupling a = 0.2 meVA, induced 
pairing potential A = 0.5 meV). 



tance for different finite values of (V x ,V y ) . It is clear 
from Fig|3] that increasing V y first suppresses the value 
of the ZBCP (at finite T — 60 mK), eventually mak- 
ing it disappear (as expected). For our chosen parame- 
ters for the system, the ZBCP essentially disappears for 
V y ~ A ~ 0.5 meV. In Fig0]we plot the actual value of 
the differential conductance at the ZBCP as a function 
of V y for V x = l meV for T = 60 and 120 mK, and it is 
clear that ZBCP would disappear for V y ~ V x , particu- 
larly at higher temperatures. We therefore predict that 
the experimentally observed ZBCP signature in Ref.[l|, 
for our estimates of the parameters of the experiment, 
should essentially completely disappear for the tilt an- 
gle 9 > 45°. An interesting notable (and experimentally 
verifiable) feature apparent in Fig. 4 is that the ZBCP is 
quite immune to a finite V y until V y becomes reasonably 
large (> 0.4 meV for our chosen parameters) when it is 
suppressed reasonably quickly. The recent experimental 
study of the ZBCP—, in fact, has studied the evolution 
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FIG. 6: Panels (a),(b) and (c) show the evolution of the 
charge density u 2 — v 2 of the split MFs generated by overlap for 
L = 4.5 pm wire with different values of applied Zeeman- field 
V x . Increasing V x increases the overlap between MFs, which 
in turn leads to a finite charge \u\ 2 — \v\ 2 7^ associated with 
the split Majorana fermion. Panel (d) shows the splitting of 
the Majorana fermion energy resulting from the same overlap 
as a function of Zeeman-field V x . The red line is calculated 
result, while the blue dashed line is a fit to the analytically 
expected overlap as discussed in the text. 



of the ZBCP as a function of transverse magnetic field 
V y , by varying the angle 8 = ta.n~ 1 (V y /V x ) in the plane 
of the wire, while holding the magnitude constant. In 
FigJS]we present our numerical results for the ZBCP as 
a function of 6 and we find oscillations of the ZBCP that 
are consistent with the experimental results in Ref.[l|. 



PROPOSED NEW EXPERIMENT 



The splitting of the MFs shown in FigJ2]depends on the 
applied Zeeman field, V x , and arises from overlap of the 
MF wave-functions localized at the two wire ends. The 
overlap of the MF wave-functions leads to Andreev states 
at non-zero energy which are no longer neutral i.e. |it| 2 — 
\v\ 2 02& in contrast to the MF mode itself, which being 
a precise zero-energy mode must be precisely a neutral 
zero-charge quasiparticle. The total ground state charge 
density in the wire can in principle be calculated using 
the expression 



</>(*)> 



E 

n:E n <0 



tanh — ^— ( I u n (x) 1 2 
2k B T Kl v ; 



k(x)| 2 ), (3) 



where the Bogoliubov quasiparticle operators with en- 
ergy E n have the form V&t = J dxu n (x)c ! (x) + v n (x)c(x). 
However, the charge density (p(x)) is influenced by the 
presence of disorder and one cannot in general separate 
the contribution of the background charge density from 
the charge density resulting from the splitting of MFs. 
The contribution of the split MFs to the charge den- 
sity can be separated from the fixed background charge 



density (which contains the contribution from impurity 
potentials) by controlling the occupancy of the split MF 
modes by coupling the ends of the superconducting sys- 
tem to an external normal fermionic lead. By chang- 
ing the chemical potential of the normal fermionic lead, 
one changes the occupancy of the split MF modes, which 
leads to a change in the charge density in the case when 
the split MF modes are charged. The change in the 
charge density as a function of chemical potential of the 
normal fermion lead, which we refer to as "non-local 
compressibility" separates out the effect of dis- 

order from the charge of the split MFs. To calculate the 
non-local compressibility, we integrate out the normal 
fermionic lead creates a self-energy for the Bogoliubov 
states of the form 



SieadM = irr z sgn(w - Mead), 



(4) 



where T is related to the transparency of the contact be- 
tween the lead and the superconductor. This leads to an 
expression for the lead chemical potential p\ ea d depen- 
dent charge density in the semiconductor 
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The non-local compressibility defined as 
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can in principle 2 ^ be measured-using single-electron- 
transistor (SET) spectroscopies in lock-in- mod o 28 i 29 . At 
low temperatures T, the compressibility singles out the 
split MF states near the chemical potential E n ~ pi ea d- 
When the value of pi ea d coincides with the energy E n , 
there would be a peak in the non-local compressibility. 
The spatial profile of this peak is related to \u\ 2 — \v\ 2 , 
which as shown in Fig |6] would show an oscillatory profile 
in space whose magnitude increases with the strength of 
the Zeeman field V x . Here one has ignored the screen- 
ing properties of the superconducting density of states. 
However, if only part of the nanowire is covered by the 
superconductor and the screening length of carriers in 
the nanowire is long compared to the oscillation wave- 
vector shown in Fig|6l the electrostatic screening effects 
of this charge density are not expected to be significant. 
FiglEId) shows the MF energy splitting resulting from 
the Zeeman splitting. This splitting (shown by the red 
curve) also oscillates as a function of the Zeeman split- 
ting V x - The splitting results from the overlap of the 
MF wave-functions^, which in the case of semiconductor 
nanowires has been shown to have the form^ ijj(L/2) = 
^e^r'CKHMV,))!^ The resulting MF overlap 
and splitting have the form splitting e~ L ^^'(a + 
bV x ) cos {k,F (V X )L + c + dV x ) where a,b,c,d are fitting- 
parameters. As seen in FigjB^d), we find the parameters 
(a, b, c, d) can be chosen to fit the V x dependence over 
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a large range. This shows that the splitting of the low- 
energy eigenstates arising from the end of the nanowircs 
shown in Figj6]is consistent with resulting from the over- 
lap of the localized end MFs. 

The above discussion and the results presented in Fig. 6 
provide an exciting possibility for a scanning SET-based 
measurement and detection of the MF mode in semi- 
conductor hybrid structures by measuring directly the 
charge associated with the split ZBCP modes in large 
Zeeman fields (where the end MFs overlap leading to 
MF-splitting), and then extrapolating this back to zero- 
splitting which should correspond to the neutral MF. 
This is by no means an easy experiment, but in prin- 
ciple, it is doable using existing SET techniques, and the 
spatially-resolved MF information coming out of such an 
experiment would be invaluable in understanding the de- 
tails not only about the existence of the MF, but also 
about its wavefunction and its spatial location. 



VI. DISCUSSION: DISORDER AND OTHER 
EFFECTS 

The theoretical work and the numerical results pre- 
sented in this work are based on a simple effective model 
which neglects many complications of the real system 
studied in the experiment. For example, although it 
is straightforward to include disorder in our model by 
adding a random one-particle potential in the tight- 
binding Hamiltonian, we have ignored providing results 
including effects of disorder because disorder effects have 
already been studied elsewhere in the literature^ and 
the details of the disorder operational in the experimen- 
tal systems are not known. The main effect of disorder 
is to suppress the topological gap (and to eventually de- 
stroy it for disorder strength comparable to the gap size) , 
and the specific experimental observations under consid- 
eration in the current work (namely, the suppression of 
ZBCP height, the Majorana splitting, the magnetic field 
effects, etc.) are unlikely to be qualitatively modified by 
disorder, particularly for the relatively clean wires stud- 
ied in Rcf. [1] . Our work also neglects any orbital coupling 
introduced by the magnetic field which could affect some 
of the experimental findings at very large fields where the 
magnetic length could be shorter than the wire width. 
Another limitation of the current theory is to approxi- 
mate the tunneling measurement simply by a single tun- 
neling barrier height U (and the finite wire length) with- 
out considering the actual details of the tunnel junction 
used in the experiments. This would necessitate a careful 
modeling of the tunnel contact which is well beyond the 
scope of the current minimal model. A full understand- 
ing of the experimental details must await a more quan- 
titative realistic modeling of the actual semiconductor- 
superconductor hybrid structure used in the laboratory, 
which may, however, be challenging since many details 
(e.g. the nature of various interfaces) are unknown. 

Two specific new results arising from our theoretical 



work are partially in agreement with observations: the 
large suppression of the ZBCP height compared with the 
pristine theoretical quantization (2e 2 /h) and the splitting 
of the ZBCP at large magnetic fields. Even for these two 
predictions, questions could be raised if alternate mech- 
anisms are operational. The strong dependence of the 
actual ZBCP on the wire length and the tunnel barrier 
that we find in our results establishes that temperature 
by itself is insufficient to explain the quantitative height 
of the ZBCP, but neither the precise effective wire length 
nor the tunnel barrier height are known, again making a 
comparison between theory and experiment problematic 
at a quantitative level. False color images can easily be 
produced by adjusting various parameters making theory 
and experiment look alike, but this may not be partic- 
ularly useful in view of the large number of unknowns 
in the problem. The Majorana splitting at high mag- 
netic fields is inevitable due to the increasing magnitude 
of the dimensionless separation between the two end Ma- 
jorana modes (as a ratio of the coherence length which 
increases monotonically with increasing field due to the 
suppression of the gap), and thus all experimental ob- 
servations of the Majorana in real systems with finite 
wire lengths must eventually reflect this splitting. The 
Majorana splitting of the ZBCP should not behave as a 
simple linear function of the magnetic field distinguish- 
ing it clearly from any trivial Zeeman splitting. A clear 
prediction of the theory is that the splitting should oscil- 
late in some manner, and such oscillations have not yet 
been seen experimentally. Whether the non-observation 
of the Majorana splitting oscillation is a genuine feature 
or is due to disorder and thermal broadening is unclear 
at this stage. Only further work can clarify this question. 



VII. CONCLUSION 

Before concluding, we point out that there are vari- 
ous resonance structures in our numerical results which 
arise from the sharpness of our confinement and trans- 
port models, which are completely non- universal and 
non-topological in nature. Such resonant structures in 
the current-voltage characteristics are well-known in ID 
systems^i and arise from various resonances in the 
transmission and reflection coefficients. Presence of dis- 
crete impurities may lead to additional non-topological 
resonant structures. These resonant structures will shift 
around with gate voltage and magnetic field with the 
ZBCP being the only universal topologically robust fea- 
ture in the data. 

We mention that although the results presented in this 
work are restricted to the one-subband strict ID limit 
(i.e. very large inter-subband gap energy) with no dis- 
order, we have carried out some representative calcula- 
tions for multisubband-occupied disordered wires finding 
qualitatively similar results, leading to our belief that 
our results and conclusions presented in this work con- 
tinue to apply qualitatively in more realistic multisub- 
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band wires in the presence of finite disorder (provided a 
TS phase can be realized in the system). While any de- 
tailed quantitative comparison between experiment and 
theory must await more realistic modelling of the actual 
SM/SC systems utilized in Ref.[l(, our current work de- 
cisively demonstrates that the suppression of the ZBCP 
well below the canonically quantized value, splitting of 
the ZBCP at high Zeeman fields, and the suppression of 
the ZBCP in the presence of a Zeeman field along the 
spin-orbit direction are all expected theoretical features 
of the SM/SC Majorana system proposed in Refs.@-(1] 
and observed in Ref.[l[. In addition, we propose a possi- 
ble new experiment using the scanning SET spectroscopy, 
which, in principle, can provide information about the 
spatial location of the MF and its wavefunction by mea- 
suring the effective charge on the Andreev bound states 
associated with the splitting of the MFs in a finite wire 
at higher values of the Zeeman field. 

Our work here is an extension and generalization of 
the standard theory for the Majorana modeSr— to finite 
wire lengths at finite temperatures. This is important 
for the comparison with the experiment 1 at a qualitative 
level since this directly tells us the extent to which the 
standard zeroth order theory captures the essential fea- 
tures of the experiment. We do not introduce any new 
ingredient except to extend the standard theory to the 
experimental finite systems at finite temperatures so that 
an honest and direct assessment is possible about the key 
issue of whether the experiment of Ref. [1] really explores 
the theoretical predictions of Refs.2-5. Such a compari- 
son with the extended standard theory is essential to es- 
tablish future directions of research in this subject. We 
do find several surprises in our results of applying the 
standard theory to finite systems corresponding to finite 
length nanowires used in the experiments. For example, 
the ZBCP never reaches its expected quantized value in 
finite length wires even at very low temperatures because 
the details of the tunneling barrier now become an im- 
portant variable. Another somewhat unexpected result, 
which should be checked in future experiments, is an os- 
cillatory splitting of the MF as a function of the effec- 
tive Zeeman splitting at high values of the magnetic field 
where the topological gap is sufficiently suppressed so 
that the coherence length is larger than the wire length. 

We conclude by stating that we have established in 
this work that the two puzzling features of the likely 
experimental observation of the Majorana modes in ID 
InSb nanowirei following earlier theoretical proposals^— 
can be explained by including finite wire length, fi- 
nite temperature, and finite tunneling barrier effects in 
the theory. We have also made specific predictions on 
how the zero-bias-conductance peak will be suppressed 
in the presence of a finite magnetic field in the trans- 
verse direction and how the Majorana-induced zero-bias- 



conductance-peak will split at high Zeeman splitting in 
a finite wire due to the overlap between the Majorana 
modes at the two ends of the wire. Our work does not 
prove that the experiment reported in Ref.[l] is indeed 
the experimental discovery of the long-awaited Majorana 
fermion, particularly since unknown effects could, in prin- 
ciple, conspire to give a weak zero bias peak which fol- 
lows for unknown reasons the phenomenology observed 
in the experiment and theoretically modeled in the cur- 
rent work. What our work does provide is compelling 
support for the claim that the observation of Ref.[l] is 
indeed consistent with the theoretically predicted exis- 
tence of the MF in semiconductor hybrid structures, and 
the apparent anomalies in the data of Ref. [1] (suppressed 
ZBCP, high-field splitting of the ZBCP, etc.) are all com- 
pletely in accord with the zeroth order theory. The ab- 
solute confirmation of the solid state MF discovery must 
await a direct demonstration of the non-Abelian nature 
of these quasiparticles through an interferometry mea- 
surement which however is unlikely to be either easy or 
quick. Meanwhile, further observations of ZBCP in other 
nanowires by different groups and various consistency 
checks between theory and experiment, as carried out 
in the current work, would go a long way in establishing 
that perhaps the elusive Majorana has finally shown up 
in a most unusual place, in a semiconductor placed on a 
superconductor in the presence of Zeeman spin splitting. 
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tum Optics Center. 

Note added: After the posting of our work, the re- 
sults of Ref. [l[ appeared on line in its published form 3 -^. 
All our results, discussion, and conclusion with respect 
to Ref.pJ apply equally well to the published results in 
Ref. [Hj]. Thus, the case in favor of the likely experi- 
mental observation 1 * 3 ^ of the possible signatures for the 
existence of the proposed" Majorana modes in SM/SC 
hybrid structures is further enhanced by our theoretical 
results. We point out, however, that at best the observa- 
tions of Refs.[Hl32| establish only the necessary conditions 
for the existence of the long-sought emergent Majorana 
modes in solid state systems. Much more work would 
be needed, including the observation of similar effects 
in other semiconductor nanowires with strong spin-orbit 
coupling (e.g. InAs) and the experimental demonstration 
of the sufficient conditions for the existence of the Ma- 
jorana modes involving the observation of the fractional 
Josephson effecl4&22 and/or the non-Abelian braiding 34 , 
before one can compellingly claim to have discovered the 
elusive Majorana quasiparticles in solid state systems. 
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